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Abstract 

The problem of controlling a single link flexible manipulator is con- 
sidered. A self-tuning adaptive control scheme is proposed which consists 
of a least squares on-line parameter identification of an equivalent linear 
model followed by a tuning of the gains of a pole placement controller 
using the parameter estimates. Since the initial parameter values for 
this model are assumed unknown, the use of arbitrarily chosen initial pa- 
rameter estimates in the adaptive controller would result in undesirable 
transient effects. Hence, the initial stage control is carried out with a 
PID controller. Once the identified parameters have converged, control 
is transferred to the adaptive controller. Naturally, the relevant issues 
in this scheme are tests for parameter convergence and minimization of 
overshoots during control switch-over. To demonstrate the effectiveness of 
the proposed scheme, simulation results are presented with an analytical 
nonlinear dynamic model of a single link flexible manipulator. 
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1 INTRODUCTION 


Automated manipulation is finding increasing use in production, 
military and space industries for performing routine, monotonous 
and hazardous tasks. The present day manipulators can perform 
with sufficiently adequate accuracy at the expense, however, of 
payload capacity and operating speed. One possible cost-effective 
solution is to build manipulators with lighter links. The lightweight 
links reduce the moment of inertia at each joint, permit the use 
of direct drive motors and have the advantages of manufacturing 
simplicity and lower cost. The next generation of manipulators 
would naturally have to be flexible. Mechanical flexiblity, however, 
generates a fairly severe problem of control of the manipulator end' 
effector motion in its work space. This is due to the inevitable 
excitation of structural vibrations and the resulting interactions 
between these vibrations and the control action which would effect 
the accuracy required of the manipulator. 

The successful implementation of flexible robots is contingent 
upon achieving acceptably uniform performance with regard to 
variations in load , task specification, reasonable speeds and the 
ability to compensate for any environmental disturbances. In con- 
trast to robots consisting of rigid links, the dynamic behavior of 
flexible manipulators is not easy to characterize, especially under 
conditions of high speed and large amplitude motion. It is not 
only the fact that this behavior is described by highly nonlinear 
differential equations but also the lack of a precise knowledge of 
this description that makes the design of an acceptable control 
system, over the total range of its operation, a formidable task. 
The dynamic effects due to changes in configuration, load and mo- 
ments of inertia, higher speed and unpredictable disturbances tend 
to degrade the performance of the flexible manipulator arm. The 
control scheme that is to be developed, therefore, has to overcome 
these dynamic effects. 

This paper attempts at a resolution of difficulties posed by this 
problem by employing a self-tuning control approach. The strategy 
here briefly consists of (i) a least squares on-line parameter identi- 
fication of an equivalent linear model, followed by (ii) a tuning of 
the controller gains by an adaptive control algorithm throughout 
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the range of the manipulator operation. Thus any changes occur- 
ing in the manipulator’s dynamic description will automatically 
be reflected in the parameter estimates and would, therefore, be 
counteracted by updating the controller gains. 

An important step prior to parameter identification is to obtain 
a valid model structure of the manipulator dynamics. This is de- 
rived by analytical modelling based on Lagrange’s equation and as- 
sumed mode shape functions from the finite element method. This 
nonlinear analytical model is used to generate the input/output 
data which, in turn, is employed in the least squares parameter 
estimation. Since, initially, the parameters are assumed to be un- 
known, the parameter estimates obtained during this initial stage 
would be unsuitable for updating the controller gains. Hence, dur- 
ing this initial stage, a simple PID stabilizing controller is used 
with the manipulator model and the parameter identification pro- 
cess is initiated. On convergence of the parameters, the control 
action is switched over to the adaptive controller. A salient fea- 
ture of the present work involves the implementation of a conver- 
gence test to minimize any undesirable transient effects following 
the switch-over. 

2 THE ADAPTIVE POLE PLACEMENT CON- 
TROLLER 

The control scheme considered here is based on adaptive pole place- 
ment. While a variety of configutations can be found in the lit- 
erature ( [1], [3], [4] ) for pole placement, the one involving a 
Luenberger observer structure (Fig. 1) as suggested by Elliot and 
Wolovich [l] is used here. This choice is based on the fact that it 
results in a closed loop system of the same order as the open-loop 
system (due to pole-zero cancellations). Also it does not add any 
undesirable zeros to the plant as might happen with the structure 
suggested in [3]. 

The adaptive pole-placement concept is briefly presented below 
in a discrete-time framework: 


3 


( 1 ) 


Let the plant to be controlled have the transfer function 

^(r 1 ) 

Aiq- 1 )' 

where q~ 1 is the backward shift operator, and 

^(g' 1 ) = l + fa,?*' 

«=i 

B(<r l ) = t>if‘ 

1=1 

so that it has the description 


( 2 ) 

(3) 

(4) 


A{q~ l )y{t) = B(q~ l )u{t). 
w here u(t) and y(t) are the input and output respectively. 

The Adaptive Pole Placement Algorithm : 

From the structure of Figure 1, one can formulate the following 
equations. 


<?(? MO = K{t,q- l )u(t) + H{t,q~ l )y(t) 


«(0 = s (0 +”(0 


where 


= Ek(t)q-* 
1=1 


«=l 


and 


Q(? 0 = i + Em"' 
1 = 1 


( 5 ) 

( 6 ) 

(7) 

( 8 ) 

(9) 
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Figure 1: Luenberger observer structure for pole placement control 


Let 

i(f,g _1 ) = 1 + (10) 

t = l 

B(t,q- l ) = 'tb i (t)q- i ( 11 ) 

i=i 

where Oi(t) and 6, (t) are the estimates of a ,• and 6,-. 

If AT(f,^* 1 ) and H{t,q~ l ) are made to satisfy the following re- 
lation 


H(t,q l )B{t,q l ) + K(t,q 1 )i(f,g l ) 

= (?(?~ l )[i(f,? _1 ) - >l<i(i,9' 1 )] 

then the resulting closed-loop transfer function becomes 

B[q^l 

A d (t,q~ l ) 


( 12 ) 


( 13 ) 
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when the identified parameters converge to the plant parameters, 
where 

Ad[t,q *) = 1 + 5"! Adj(t)q~' (14) 

i=i 

With this structure, however, the plant cannot be made to track 
a step input signal. In order to equip this structure with such a 
tracking facility, unity feedback is applied and an integrator is 
inserted in the forward path. This can be formulated as 


w ( t ) = >s(0 • c + w(t - 1) (15) 

s(t) = v(t) - y(f) (16) 

Also equation ( 6 ) should be modified to : 


u{t) = g[t) + w(t) 


(17) 


The desired denominator A d (t,g _1 ) and the scalar gain c can 
be determined from the desired closed-loop denominator £?(y -1 ) 
by the following equations. 


1 

-1 

0 

0 

0 


0 0 0 • • ■ \ 

i o o ■ • • S 2 

-11 0 ••• 6s 

0 0 -1 1 6 n 

0 0 0 -10 


( Adl ) 


( d\ + 1 ^ 

Ad2 


d 2 

Ad 3 

— 

dz 

A-dn 


i 

\ c ) 


[ dn+l 7 


where 


n+l 

D(q~ 1 ) = 1 + E *9-* 
1 = 1 


(18) 


(19) 


Since we can obtain Ad{t,q~ l ) from equation (18) , the H{t y q~ l ) 
and .RT(£,g -1 ) can be obtained from equation (12). 

The block diagram of this scheme is shown in figure 2. 

Note that the step input tracking facility is achieved by increas- 
ing the order of the overall system to only n+l. 
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Figure 2: Adaptive Pole Placement Control Scheme Using Luenberger Observer 
Structure with Integral Action 


The Least Squares Identification Algorithm: 

The estimates a* and 6 t used in the control scheme are obtained 
by a least squares parameter identification algorithm [4] as follows: 


0{t) = 6{t — l) + 


l + $(t-l) r P(t-2))3(t-l) ’ K } 


t > 1 


P{t - 1) = P{t - 2) - 


P{t - 2 )$(t - 1 )$(* - 1 ) T P{t - 2) 
1 + $(t - 1 ) T P(t - 2)$(t - 1) 


( 21 ) 


with 0(0) given and P(— 1) is any positive definite matrix, 
where 


6(t) = (-di(t), -a 2 (t),- • • ,-a n (t),Si(t),S 2 (t), - •• ,6 n (t)] T 
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is the current parameter estimate vector, and 


$(t-l) _ [y(f _ l), y[t — 2), • • • ,y{t - n),u(t - l),u(f - 2),---,u(t - n)] r 


3 SWITCH-OVER FROM PID TO ADAPTIVE 
CONTROLLER 

A critical question in the present control scheme is to determine 
an appropriate time to switch from the initial PID controller to 
the adaptive controller. The simplest way is to wait till the pa- 
rameter estimates resulting from the identification algorithm have 
converged to their true values. The following criterion provides a 
check on such a convergence. 

The Convergence Criterion: 


Assume || 0(0) || 2 < M . 
If 


then 

where 


l)] ^ 


^m»n[P( l)] 
M 


( 22 ) 


ii m ii 2 < 


(23) 


Oft) = $(t) - 0 

0(t) : identified parameter vector at t im e t 
6 : actual parameter vector = 

[ — a l> ~ a 2 * ' • » 

e : the error tolerence for the convergence test 
of the identified parameters. 

^max[.P(t - 1)] : maximum eigen value of P(t - 1) . 
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Proof: 

From [4] (p. 61), one can get the following inequality 


A m ,n(p(* - inn m ii 2 < A mai (p(-i)- i ni no) ii 2 

which implies 

11 5 11 * (0 > »* 

Using (22) and (25), it follows that 


m n* < ^ 


< € 


(24) 

(25) 

(26) 


Thus by computing A max [P{t — 1)], one can test the convergence 
of the parameter estimates. 

Switching Logic: 

Once the identified parameters have converged to their true val- 
ues and the system step response has reached steady state, control 
action is switched over from the PID to the adaptive controller. 
This is probably the simplest manner to implement the switch- 
over without causing any undesirable transients. 

An alternative switching logic is proposed here which does not 
require the step response to reach steady state. However, this 
logic is limited only to those systems that satisfy the conditions 
for one-step-ahead control [4]. 

Assume 

(i) the plant to be linear time invariant, 

(ii) the switching instant to correspond to t = 0, and 

(iii) the desired output trajectory after switching to be the 
same as the one that would have been obtained, had the adap- 
tive controller been applied to the plant starting at rest from that 
position(y 0 ), where 

y a = [y( — l),y(— 2), • • • ,y( — n),u(— l),u(— 2),- • • ,u(— n)] r 0 (27) 
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is the output at the switching instant. 

In order to satisfy the last assumption, a correction input u e (£) 
is needed to compensate for the terminal conditions resulting from 
the PID controller. Thus the plant input would be 

u(«) = u d {t) + u c (t) (28) 

where u d (t) is the input generated by the pole placement algorithm. 

The plant output can be expressed as 

y(0 = Vd{t) +y e {t) (29) 

where 

I Vd(t) = $/(t)0 ; y<<(0) = 0 (30) 

y«(0 = (0^ iyc(o) = y 0 (3i) 

where the subscripts denote the correspondence of the two compo- 
nents. 

From assumption (iii), j/c(t) = y 0 for t > 0, and the compensat- 
ing input u c (f) is obtained using (31) as: 

tt «(0 = ^■{yo-[yc(O>yc(<-i),---,y< : (<-»+i),o,u e (i-i),...,u c (t-n+i)] T 0} 

(32) 

With a proper choice of the sampling interval, the flexible ma- 
nipulator discrete model is found to meet the requirements of one- 
step-ahead control. However, this approach is found to be suitable 
only in those situations where the deflections are small, and is not 
used in the simulation here. 

The complete control block diagram is shown in Fig. 3 . 

4 SIMULATION RESULTS 

The dynamic analytical model of the single link flexible manipula- 
tor is described by [5]: 
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lire 3: Complete Block Diagram of PID-Adaptive Pole Placement Controller 
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{J + Eili m it r?) 
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m 2 p 

m N0 


where 


m l p ^ 


0] 


u ~ 2 E.^i m* ;r,r, \ 

0 0 

N 

0 

H 

s 0 


r~i 

f 2 

— 

(0 2 ~ul)ri 
{0 2 - w|)r 2 

: ••• i 0 

0 0 ... m NN J 


V. r~N J 


l (0 2 — wjf)r N J 


(33) 

y = tan — j-!- + 0 (34) 


N 

*»(*.<) ®£&(*) r <(<) (35) 

«=i 


u : input torque to the beam 
y : tip position 
/ : the length of the beam 
0 : the hub angle 

4>i{x)\ the mode shape functions of the beam 
r,(t) : the generalized coordinates 

The desired closed-loop denominator for the adaptive pole place- 
ment controller is chosen as 


D{q l ) = 1 + d x q 1 + d 2 q 2 (36) 

where 


d\ = -2 e~ iujT cos uTyJl - £ 2 (37) 

d 2 = e~ 2€wr (38) 

where T is the sampling period in seconds. 

For computer simulation, the following numerical values are 
used: n = 4, N = 2, u = 5, £ = 1, T = 0.1, P(-l) = 10 8 / 8x8 , 
M = 10 and e = 0.7. The switching from PID to the adaptive 
controller occurs at t = 4 secs. 

The results are shown in Fig. 4. 
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Figure 4: Simulation of Combined PID and Adaptive Pole Place- 
ment Control: 

(A) reference input v(t) and plant output y(t), 

(B) plant input u(t), 

(C) convergence of identified coefficients of denominator, 

(D) convergence of identified coefficients of numerator. 
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5 CONCLUSION 


A simulation based study for the adaptive control of a single link 
flexible manipulator has been considered. Such a control approach 
is of practical importance since the dynamic characteristics of the 
manipulator change considerably especially while picking up or re- 
leasing payloads. In such cases unless the control gains are suitably 
updated, the performance would be poor. 

Since the adaptive control scheme depends on the parameter 
estimates from an on-line identification algorithm, the initial con- 
trol action is carried out with a PID controller during which the 
identification process is initiated. On convergence of the parame- 
ter estimates, control is smoothly transferred to the adaptive con- 
troller. A criterion for testing the convergence has been presented. 
The simulation results amply demonstrate the effectiveness of the 
proposed scheme. 

Experimental verification of the control scheme on a laboratory 
test set-up is presently in progress. 
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